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ABSTRACT 

We present the results of a search for cool white dwarfs in the United King- 
dom InfraRed Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS) Large 
Area Survey (LAS). The UKIDSS LAS photometry was paired with the Sloan 
Digital Sky Survey (SDSS) to identify cool hydrogen-rich white dwarf candidates 
by their neutral optical colors and blue near-infrared colors, as well as faint Re- 
duced Proper Motion magnitudes. Optical spectroscopy was obtained at Gemini 
Observatory, and showed the majority of the candidates to be newly identified 
cool degenerates, with a small number of G- to K-type (sub) dwarf contaminants. 
Our initial search of 280 deg 2 of sky resulted in seven new white dwarfs with 
effective temperature T e g « 6000 K. The current followup of 1400 deg 2 of sky 
has produced thirteen new white dwarfs. Model fits to the photometry show 
that seven of the newly identified white dwarfs have 4120 K < T c s < 4480 K, 
and cooling ages between 7.3 Gyr and 8.7 Gyr; they have 40 kms -1 < v tan < 85 
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kms -1 and are likely to be thick disk 10-11 Gyr-old objects. The other half of the 
sample has 4610 K < T c g < 5260 K, cooling ages between 4.3 Gyr and 6.9 Gyr, 
and 60 kms" 1 < v tan < 100 kms -1 . These are either thin disk remnants with 
unusually high velocities, or lower-mass remnants of thick disk or halo late-F or 
G stars. 

Subject headings: white dwarfs — techniques: photometric — techniques: spec- 
troscopic — Infrared: Stars — surveys 

1. Introduction 

White dwarfs are the end stage of stellar evolution for the vast majority of stars — all 
stars with initial mass between around O.O7M and 8M end their lives as cooling white 
dwarfs. The cooling rate of these degenerate remnants slows as their temperature drops, such 
that very old white dwarfs are still visible. The coolest white dwarfs can constrain the age 
of the Galactic disk, or even of the halo (e.g. Winget et al. 1987). While other methods for 
dating the Galactic components exist — asteroseismology, gyrochronology, isochrone fitting, 
isotope decay, magnetic activity, kinematics and metallicity (e.g. Ulrich 1986, Barnes 2007, 
Chaboyer et al. 1996, Frebel et al. 2007, West et al. 2008, Nordstrom et al. 2004, see also the 
review Soderblom 2010) — very cool white dwarfs are unambiguously old and their cooling 
rates are quite well understood, enabling their use as accurate chronometers (e.g. Iben & 
Tutukov 1984; Fontaine, Brassard & Bergeron 2001). Most white dwarfs consist of a C/O 
core with an outer envelope composed of helium and/or hydrogen, with occasional traces 
of metals. The mass and composition of both the core and the atmosphere are important 
in determining the age of a white dwarf. The mass and core composition determine the 
thermal content, while the composition of the atmosphere provides the insulating material 
which controls the rate of cooling. Bergeron et al. (2001) use atmospheric and evolutionary 
models to analyse a sample of white dwarfs with measured trigonometric parallaxes to show 
that the coolest of these white dwarfs, with T ef j ~4000-4500 K, are 9-10 Gyr old if they 
have a thick hydrogen atmosphere, and 8-9 Gyr old if they have a helium-rich atmosphere. 
These ages are consistent with the age of the local Galactic disk (e.g. Leggett et al. 1998). 

Several groups are trying to find old white dwarfs to confirm the age of the disk and to 
investigate the ages of the older Galaxy components. Candidate white dwarfs are typically 
identified as high proper motion objects, or by their colors or spectra, or by a combination 
of their kinematic and photometric properties. Using kinematic data alone is problemati- 
cal. Oppenheimer et al. (2001) identified a sample of high- velocity white dwarfs which was 
inferred to be a halo population by their kinematics. However Reid et al. (2001) suggest 
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that the majority of this sample has kinematics consistent with thick disk membership, and 
analysis of the sample by Bergeron et al. (2005) found that the white dwarfs were relatively 
warm, implying relatively short cooling ages. The short cooling age does not necessarily ex- 
clude the possibility that these white dwarfs are old — they may have evolved from low-mass 
stars with long main-sequence lifetimes. 

Omitting infrared data in studies of candidate cool white dwarfs can also lead to large 
uncertainties in the derived T e g and hence age. This is because the complete optical to 
infrared spectral energy distribution is required to determine the atmospheric composition 
and T e ff. Hydrogen-rich white dwarfs cooler than about 5000 K (depending on instrument 
resolution) are featureless, however the infrared region demonstrates the presence of hydrogen 
due to the opacity of pressure- induced H 2 absorption (e.g. Borysow 2002). Kilic et al. 
(2010a) perform a detailed analysis of cool white dwarfs using optical spectroscopy and 
infrared photometry. They determine significantly different values of T e s from those found 
using optical data only; several white dwarfs thought to have T e g ~3500 K in fact have 
T e ff ~4500 K, and the cooling age decreased from >10 Gyr to <8 Gyr. 

Very cool white dwarfs are unambiguously old, as their total age is dominated by the 
large, post-main sequence, cooling time. There are about 20 white dwarfs known with 
T c ff < 4000 K. These have been found as high proper-motion objects in photographic sky 
surveys (Hambly et al. 1997, Harris et al. 1999, Ibata et al. 2000, Oppenheimer et al. 2001, 
Rowell et al. 2008, Ruiz & Bergeron 2001, Scholz et al. 2002) or in the Sloan Digital Sky 
Survey (SDSS, York et al. 2000) by inspection of spectra or by using the Reduced Proper 
Motion (RPM; Gates et al. 2004 , Hall et al. 2008, Harris et al. 2001, Harris et al. 2008, 
Kilic et al. 2006, Kilic et al. 2010b). The RPM is defined as: 

H = mag + (5 x log(PM) + 5) 

where PM is the proper motion in " yr -1 , and the RPM acts as a proxy for absolute 
magnitude for a sample with similar kinematics. Faint white dwarfs with significant proper 
motion can be separated from stars and subdwarfs using an RPM diagram, such as that 
shown in Figure 1. 

Our group has paired the UKIRT Infrared Deep Sky Survey (UKIDSS, Lawrence et 
al. 2007) with the SDSS, to provide a full optical through near-infrared spectral energy 
distribution, and to select candidate cool white dwarfs using the RPM and various colors. 
The proper motions are derived using the target coordinates and the epochs of the SDSS and 
UKIDSS images. Combining these databases allows us to go about a magnitude deeper in 
the g-band than other searches, to g « 21, although we are also limited to r « 21 by the time 
required to obtain a spectrum for classification. In Lodieu et al. (2009b) we published the 
results of a search of 280 deg 2 of the Large Area Survey (LAS) component of the UKIDSS 
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Data Release (DR) 2. The search produced seven new white dwarfs with T e g m 6000 K, 
which were confirmed with optical spectroscopy obtained at Gemini Observatory. Here we 
report the identification and confirmation of white dwarfs selected from 1400 deg 2 of DR 6 
of the LAS. The larger area and better understanding of the photometry has allowed us to 
discover significantly cooler and older white dwarfs compared to the Lodieu et al. sample. 



2. Sample Selection 

2.1. The LAS and SDSS Databases 

UKIDSS consists of five survey components, one of which is the LAS (Lawrence et al. 
2007). The LAS is the sub-survey most likely to contain faint and rare sources of the local 
Galaxy, such as the cool white dwarfs and brown dwarfs. The LAS aims to survey 4000 
deg 2 in YJHK (Hewett et al. 2006) with a second epoch at J, to reach J ~ 20. The 5a 
photometric depths are Y = 20.2, J = 19.6, H = 18.8 and K = 18.2 (Dye et al. 2006). The 
data and catalogs generated by the automatic pipeline processing can be retrieved through 
the WFCAM Science Archive (WSA; Hambly et al. 2008). All data are pipeline-processed 
by the Cambridge Astronomical Survey Unit (CASU) following a standard procedure for 
infrared images (Irwin et al. 2004). An extensive description of each step involved in the 
processing of the WFCAM data is available on the CASU webpag^j]. 

The work presented here involves spectroscopic follow-up during Gemini Observatory 
semesters 2008B and 2010A (see §3). Our selection of candidates evolved over time, as 
LAS Data Releases became available, and as we obtained spectra for the candidates. The 
Observatory was flexible in allowing us to replace targets listed in the proposals with better 
candidates as they were found in new Data Releases. The LAS Releases used were DR3 
(December 2007, 809 deg 2 ), DR4 (July 2008, 984 deg 2 ), DR5 (April 2009, 1270 deg 2 ) and 
DR6 (October 2009, 1434 deg 2 ). 

We have cross-correlated the LAS with the SDSS database, using Structured Query Lan- 
guage (SQL) and the WSA. SDSS Data Releases 6 and 7 were used, as available (Adelman- 
McCarthy et al. 2008, Abazajian et al. 2009). The area surveyed by the LAS was designed 
to overlap with the SDSS, divided up into three blocks (Dye et al. 2006, Lawrence et al. 
2007, Warren et al. 2007). The equatorial block with Right Ascension 22 to 04 hours and 
Declination between —1.5 and +1.5 degrees overlaps SDSS stripes 9 to 16. The southern 
block covers 8 to 16 hours and (approximately) —3 to +15 degrees, and includes SDSS stripe 
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82. The northern block covers 8 to 17 hours and (approximately) +30 to +50 degrees, and 
provides an overlap with SDSS stripes 26 to 33. 

2.2. Astrometry and Photometry 

Sources were matched by requiring the presence of a "primary" SDSS source within 5" of 
the LAS coordinates (increasing the search radius led to erroneous pairings). The WFCAM 
astrometry is tied to the 2MASS point source catalog (Skrutskie et al. 2006) and has a 
systematic accuracy of < O'.'l rms (Dye et al. 2006). Lodieu et al. (2009) plot the rms of 
the difference between the coordinates of all point sources in the LAS DR2 and SDSS DR5 
as a function of J magnitude; for the sources in this sample with 17.5 < J < 19.2 the implied 
astrometric uncertainty is 0'/025 to 0'.'035 for the brighter to fainter sources, respectively. 

We have restricted our queries of the LAS database to detections classified as point 
sources, and imposed color criteria and lower limits on the g-band RPM, H g . The queries 
return coordinates, photometry and errors from both surveys, as well as the proper motion, 
computed from the difference between the LAS J-band coordinates and the SDSS z-band 
coordinates. For the sample presented here, the SDSS epoch ranges from 1999 to 2006, and 
the LAS from 2005 to 2008. The average time period between the LAS and SDSS astrometry 
is 4.5 ± 2.0 years. Given the astrometric uncertainty, the uncertainty in proper motion is < 
14 mas. For one source, ULAS J1323+12, the LAS and SDSS epochs are separated by only 
0.89 yr and the proper motion uncertainty is 30 mas. 

Table 1 gives astrometric information for the candidate white dwarfs studied here, and 
Table 2 gives the SDSS ugriz and the LAS YJHK photometry. Note that the photometry 
is taken from the 8th Data Releases for both the LAS and SDSS, and not the releases used 
for our initial source identification. The proper motion in Table 1 is similarly an updated 
value, and has been calculated from the coordinates and epochs given in the DR8 releases. 

With the exception of the close binary ULAS J1234+06, Table 1 also lists proper motions 
derived from alternative catalog matching. For 10 of the 16 sources, proper motion derived 
by matching the USNO-B catalog (Monet et al. 2003) to the SDSS catalog is available 
via DR8 of the SDSS. Of these ten, there is good agreement for six objects. For two the 
differences are at the ~ 2a level (ULAS JO 12 1-00 and ULAS J 1142+00). For the remaining 
two objects the discrepancies are large, and we rederived the proper motions using the 
digitised sky images, as well as the USNO-B, SDSS and LAS astrometry, in order to obtain 
the motion over a large timeline. For ULAS J1323+12, which has a small SDSS-LAS epoch 
difference, the derived motion is close to the USNO-SDSS value, and significantly smaller 
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than the SDSS-LAS value. For ULAS J2246-00, the derived value is close to the SDSS-LAS 
value but significantly different from the USNO-SDSS value in direction. For the six objects 
without USNO-SDSS proper motions, three show good agreement between the SDSS-LAS 
proper motions and new measurements derived from the imaging and catalog data listed 
above. Three objects differ by ~ 2a (ULAS J0840+05, ULAS J1345+15, ULAS J1454-01). 
Table 1 gives the RPM values derived from both proper motion measurement sets. 

2.3. Color and RPM Selection 

Figures 1 and 2 illustrate the selection criteria imposed to arrive at the candidate list of 
Table 1. For all queries, we selected 14 < J < 19.6, to avoid saturated sources, and to ensure 
a good (5a) J-band detection, in the LAS. We also restricted the targets to tab < 20.7, so 
that the spectroscopic observations were of reasonable duration (see §3.1). To select for faint 
high-proper motion sources, we restricted sources to H g > 20.5 (see Figure 1). This last 
selection greatly reduced the number of objects returned from a query. Figure 3 shows the 
red to far-red colors of various samples, for completeness. 

Most sources, 9 of the 17 listed in Table 1, were found in the search identified as region 
A in Figures 1 and 2, and in Table 1. The cuts imposed, in addition to those listed above, 
were: H < 18.9, 0.8 < g - r < 1.6, 0.2 < r - % < 0.6, 0.6 < % - J < 1.4, J - H < 0.2. 
This selection is designed to produce hydrogen-rich white dwarfs with T eS ~ 4000 K (Figure 
2), with a good detection in the if-band. Early on in this work an additional source on the 
blue side of this cut with g — r = 0.63 was followed up; later we prioritized redder (in g — r) 
sources in order to find white dwarfs cooler than 5000 K. This bluer source is identified as 
found in region Ab in Table 1. An additional 6 sources were found in search B, defined as: 
0.2 < g-r < 1.2, -0.6 < r-i < +0.6, J-H < -0.1, H-K < -0.1. This search is designed 
to find bluer sources where H 2 opacity is impacting the red as well as the near-infrared, due 
to either lower T eS , or a higher-pressure mixed H-He atmosphere (Figure 2). Finally, one 
more source was found by selecting for extreme sources in the RPM diagram (Figure 1): 
H g > 22.5, -0.5 <g-i< 1.0, 0.0 < r - i < 0.6, H < 18.8, J-H < 0.1, H -K < 0. This 
source is identified as found in region C in Table 1 (later the RPM value was found to be 
erroneous) . 

Our spectroscopic followup is complete for the pairing of the 1270 deg 2 of the LAS DR5 
and the SDSS DR7 with the following selection: H g > 20.5, r < 20.7, 14 < J < 19.6, H < 
18.9, 0.9 < g-r < 1.6, 0.2 < r-i < 0.6, 0.6 <i — J < 1.4, J-H+ ( ^err( J) 2 + err(H) 2 ) < 
0.2. This search of the 1270 deg 2 of LAS DR5 produced six new white dwarfs with 4120 K 
< T c s < 4380 K (see §4), and recovered a known white dwarf with T e s = 4390 K (SDSS 
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J115814.52+000458.3, Kilic et al. 2010a). Kilic et al. (2010a) perform a detailed analysis 
of 113 white dwarfs with t> tan > 20 kms" 1 found using an RPM search of the 5282 deg 2 
footprint of the SDSS DR3 (Harris et al. 2006). We can compare our results to those of 
Kilic et al., neglecting the additional 13 white dwarfs in the Kilic et al. sample that were 
found through a search of SDSS spectral data, in order to compare similarly determined 
samples. Kilic et al. find 6 white dwarfs with 4100 K < T cS < 4400 K out to a distance of 
about 55 pc. The implied space density of 4200 K white dwarfs, with cooling age 8 Gyr 
(§4), is 2 x 10~ 5 pc~ 3 . The r-, J- and if-band cuts imposed for our LAS search indicate 
that we should be sensitive to 4200 K white dwarfs at distances of 10 - 95 pc, based on the 
models. The H g selection combined with a faint limit of g ~ 21 implies that our lower limit 
on proper- motion is 0'/08 yr _1 . Our upper limit is determined by the 5" matching radius 
and the difference between the SDSS and LAS epochs; the typical difference is four years 
implying an upper limit to the proper motion of ~ 1'.'3 yr _1 . Excluding ULAS J1142+00 
and ULAS J1234+06AB for which erroneously high proper motions were initially derived, 
the sample has a range in proper motion of / /09 to 0'/3 yr _1 . The density of 4200 K white 
dwarfs implied by the Kilic et al. sample suggests that we would expect to find around seven 
4200 K white dwarfs. However two of the six Kilic et al. white dwarfs have H g < 20.5, 
suggesting that our selection would find around five 4200 K white dwarfs, very close to the 
number found. 

We note in passing that we have applied this last selection to the 48 deg 2 overlap area of 
the UKIDSS Galactic Cluster Survey (GCS) DR8 and SDSS DR7, and found no candidates. 
The GCS will survey ten open star clusters in the ZYJHK filters, to a similar depth in 
YJHK as the LAS. It has been successfully used to find field brown dwarfs (Lodieu et al. 
2009a), as well as to study open clusters (Lodieu et al. 2007), but a larger area in common 
with the SDSS is required before cool white dwarfs will be found. 

3. Observations 

3.1. Optical Spectroscopy 

We obtained optical spectra of the candidate cool white dwarfs listed in Table 1, selected 
as described in the previous Section, excluding SDSS J1247+06 for which Kilic et al. 2010a 
present an optical spectrum (a spectrum is also available in the SDSS database). The 
Gemini Multi-Object Spectrographs (GMOS, Hook et al. 2004) at both Gemini North and 
South were used, through queue time granted under programs GS-2008B-Q-35, GN-2008B- 
Q-lll and GS-2010A-Q-58. For all observations the R400 grating was used with the GG455 
blocking filter. The central wavelength was 680 nm, with wavelength coverage of 460 - 890 
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nm. Detector fringing affected the spectra longwards of about 800 nm. We did not correct 
for the fringing as this wavelength region is not necessary for target classification. The 0'.'75 
slit was used with 2x2 binning, and the resulting resolution was R « 1280 or 6 A. An 
observing log is given in Table 3. 

Flatfielding and wavelength calibration were achieved using lamps in the on-telescope 
calibration unit. Spectrophotometric DA white dwarf standards were used to determine the 
instrument response curve, and flux calibrate the spectra; for program GS-2008B-Q-35 LTT 
7987 was used, for GN-2008B-Q-111 G191-B2B was used, and for GS-2010A-Q-58 LTT 3218 
was used. The data were reduced using routines supplied in the IRAF Gemini package. 

Figure 4 shows the GMOS spectra of the three non-white dwarf objects in our sample, as 
well as G- and K-type stellar spectra taken from the spectral atlas of Le Borgne et al. (2003) 
for reference, for a wavelength range that contains the features useful for classification. The 
narrow metal lines seen in the spectra exclude the possibility that these are white dwarfs. 
The three objects (ULAS Jl 142+00, ULAS J1234+06N and ULAS J2246-00) are early-G 
to early-K dwarfs or subdwarfs. For two of these, ULAS J 1142+00 and ULAS J1234+06N, 
we initially derived too high a proper motion, the revised smaller H g value places them near 
the subdwarf region of the RPM (Figure 1). Our true contamination is therefore small, one 
in fifteen, or 7%. 

Figure 5 shows the spectra for the thirteen newly confirmed white dwarfs. Two of these 
show weak pressure-broadened H a (ULAS J1323+12 and ULAS J2331+15), and eleven are 
featureless. 

3.2. Infrared Photometry 

When carrying out our analysis by fitting models to the observed photometry (§4), we 
found several instances where the infrared photometry appeared to be in error, in particular 
for the if-band. To trace the source of the discrepancies, we obtained repeat photometry 
for a subset of the sample. 

H and K photometry for five of our targets was obtained using the Wide Field Cam- 
era (WFCAM, Casali et al. 2007) on UKIRT, through the UKIRT Service program, via 
program USERV1876. Single-pointing WFCAM "paw-print" observations were defined to 
place the target in one of the four cameras: camera number 2. Exposure times of 10 s were 
used, with an 8 position telescope dither pattern, and a 4 position microstep pattern. These 
32 exposures were repeated three times, for a total on-source time of 16 minutes, in each 
filter. An observing log is given in Table 4, together with the derived photometry. The 
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"flat-file" FITS images and catalogs produced by the CASU pipeline were accessed through 
the WSA (§2). The 2"-diameter aperture photometry, provided in counts, were converted to 
magnitudes using the airmass and zeropoints provided for each of the three reduced groups 
of observations for each filter and each object. We derived the weighted mean (using the 
uncertainties provided in the catalog) and the uncertainty in the mean for these measure- 
ments. Subsequently, the merged catalog data was released for the three objects observed 
on 2010 June 13, and the photometry from this catalog, database USERV1876v20101104, is 
given in the Table for ULAS J1206+03, ULAS J1351+12 and ULAS J1404+13. The merged 
catalog photometry agreed with the flatfile photometry to within 0.02 magnitudes for the 
two brighter sources, within 3% at H for the fainter ULAS J1404+13, and within 12% at 
K for ULAS J1404+13, which has K 19.7. These differences are well within the quoted 
uncertainties. 

Additional H, or H and K, photometry was obtained for five more white dwarfs using 
NIRI (Hodapp et al. 2003) on Gemini North, via program GN-2011A-Q-59. Exposure 
times of 15 s or 30 s were used, with a 5 or 9 position telescope dither pattern. The total 
integration time is given in Table 4, together with the derived photometry. The data were 
reduced using routines supplied in the IRAF Gemini package. UKIRT Faint Standards were 
used for calibration (FS 16, FS 33, FS 126, FS 136; Leggett et al. 2006). 

Comparing Table 4 to the LAS photometry in Table 2, shows that typically, for the 
sources fainter than H m 18, the LAS if-band magnitude is too faint by about 2a. In 
selecting faint sources near the LAS detection limit, with blue J—H colors, objects have been 
scattered into the sample with H magnitudes too faint by twice the estimated uncertainty. 
A similar effect is seen in searches for faint brown dwarfs that are also blue in J — H (B. 
Burningham, private communication, 2010). 

4. Analysis 

4.1. Model Fitting 

The pure-hydrogen, pure-helium and mixed composition model atmospheres used in 
this analysis are described in detail in Kilic et al. (2010a). These models are in local 
thermodynamic equilibrium, they allow energy transport by convection, and they can be 
calculated with arbitrary amounts of hydrogen and helium. For the sample considered here, 
the mixed composition atmospheres had ratios of He to H by number ranging from low values 
of 10~ 2 to high values of 10 10 . The method used to fit the photometric data is described at 
length in Bergeron et al. (2001). Briefly, we convert the magnitudes into observed fluxes 
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using the method of Holberg & Bergeron (2006) and the appropriate filters. We use the 
SDSS to AB system corrections for u, i and z given in Eisenstein et al. (2006, Section 2). 
Then we fit the resulting energy distributions with those derived from model atmosphere 
calculations, using a nonlinear least-squares method. Only T e g and the solid angle ir(R/D) 2 , 
where R is the radius of the star and D is its distance from Earth, are considered free 
parameters. Since no parallax measurements are available, we assume a surface gravity of 
logg = 8 which determines the value of R for a given value of T cff (Bergeron et al. 2001). 
White dwarfs have been shown to have a very strongly peaked mass and surface gravity 
distribution (e.g. Bergeron et al. 1992; Liebert, Bergeron & Holberg 2005; Kepler et al. 
2007). DA white dwarfs have a mean mass of 0.6 ±0.1 M Q while DBs are slightly more 
massive with 0.7 ± 0.1 M Q ; these ranges infer a likely range in gravity for our sample of 
7.7 < log g < 8.3. The photometric variance uncertainties in T e g and the solid angle are 
obtained directly from the covariance matrix of the fit. 

Since our models do not include the red wing opacity from Ly a calculated by Kowalski 
& Saumon (2006), we neglect here the u bandpass in our fitting procedure as well as the 
g bandpass for white dwarfs cooler than 4600 K in hydrogen-rich solutions (as the missing 
Ly a opacity has a larger impact at lower temperatures). The Ly a opacity affects a wave- 
length region where there is very little flux; hence the atmospheric structure is not affected 
significantly by a change of the opacity in the ultraviolet, although the predicted blue colors 
are. 

As described in §3.2, we have found that in selecting faint sources near the LAS detection 
limit, with blue J — H colors, objects have been scattered into the sample with H magnitudes 
too faint by twice the estimated uncertainty. In our fits to the data we use the NIRI and 
WFCAM photometry where available (Table 4), and where LAS photometry is used we 
exclude data with uncertainties >0.15 magnitudes. The fits use the DR8 releases of both 
the SDSS and LAS photometry, the most recent at the time of writing. 

We use the energy distributions together with the optical spectra at H a to constrain 
the surface composition. Only two of the white dwarfs show H a (Figure 5), however this 
does not necessarily imply that the remaining objects are helium-rich, as H a absorption is 
not seen in hydrogen-rich white dwarfs cooler than ~ 5000 K. Figures 6 though 11 show the 
model fits to the observational data, assuming \ogg = 8.0. 

The energy distributions are somewhat sensitive to log g in the temperature range con- 
sidered here. A variation of ±0.3 dex in \ogg yields differences in effective temperature of 
between 40 K and 80 K. The uncertainty due to the photometric variations is between 40 K 
and 150 K. Another error estimate is provided by fitting the observed H a line profile for 
the two DAs: the derived T e g- differs from that derived from the energy distribution by 100 
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- 200 K (for ULAS J1323+12 5400 K cf. 5260 K, for ULAS J2331 + 15 5200 K cf. 4970 K). 
This is consistent with our other error estimates. 



4.2. Derived Properties 

Table 5 lists the derived atmospheric properties of the 13 new white dwarfs in our sample, 
together with the SDSS white dwarf recovered here. Using the composition and temperature, 
and assuming that these stars have the canonical white dwarf mass of 0.6 M , we can use 
the synthetic colors of Holberg & Bergeron (2006, an extension of Bergeron, Wesemael & 
Beauchamp 1995) and the evolutionary sequences of Fontaine, Brassard & Bergeron (2001) 
to derive both a cooling age and distance, and hence tangential velocity. These values are 
also given in Table 5. The derived tangential velocities, H g and g — % are consistent with 
the 40-150 kms -1 region indicated in Figure 1. The uncertainty in the implied cooling age, 
distance and velocity are primarily due to the uncertainty in gravity (or mass). The sense 
of the gravity effect is that more massive white dwarfs will have a longer cooling age and be 
closer and slower, and vice versa. 

Four of the thirteen newly identified white dwarfs (~ 30%) are best fit with pure- 
hydrogen atmospheres (Figures 6 and 7; ULAS J0826-00, ULAS J1323+12, ULAS J1345+15, 
ULAS J2331+15), and four (~ 30%) are best fit with pure-helium atmospheres (Figure 8; 
ULAS J1006+09, ULAS J1206+03, ULAS J1351 + 12, ULAS J1454-01). An additional white 
dwarf (ULAS J1320+08) is fit with either pure-helium or a helium-rich atmosphere, where 
the latter fit, with H/He = 10" 44 , is superior (Figure 9). Two sources (ULAS J0840+05, 
ULAS J1404+13) have mixed atmospheres with large flux deficits in the infared, due to 
pressure-induced H 2 absorption (Figure 10). As described in Kilic et al. (2010a), there are 
generally two solutions to such objects. The H 2 opacity peaks around H/He = 10~ 2 and 
there is usually a good solution both above and below this peak with slightly different tem- 
peratures. This is demonstrated by the fits shown in Figure 10, where ULAS J0840+05 is 
fit with H/He = 10~ 4 - 5 or 10~ a6 , and ULAS J1404+13 is fit with H/He = 10~ 2 - 8 or lO" 21 . 
Thus three white dwarfs (~ 25%) have mixed atmospheres, and are similar to the mixed 
composition white dwarfs found by Kilic et al. (2010a, their Figure 14). There are seven 
white dwarfs (~ 55%) that have pure-helium or helium-dominated atmospheres. Finally, 
there are two white dwarfs (~ 15% of the sample) for which we cannot constrain the atmo- 
spheric composition: ULAS J0121-00 and ULAS J1436+05 (Figure 11). ULAS J1436+05 
is particularly puzzling — it appears to be much too bright at zY J. A binary solution does 
not seem likely given the flux distribution, and the photometric uncertainties are small. 
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4.3. Discussion 

Typically, studies of samples of cool white dwarfs find a much larger fraction of pure- 
hydrogen atmospheres: 50-60% (Bergeron et al. 2001, Kilic et al. 2010a). However Kilic et 
al. (2010a) find that ~ 60% of their sample of white dwarfs (25 of 40) in the temperature 
range 4500-5000 K have pure-helium atmospheres while only ~ 20% have pure-hydrogen, 
consistent with our findings. As Kilic et al. state, further work is required to understand 
if the observed overabundance of helium-rich atmosphere white dwarfs at this temperature 
range is real, or if Ly a or H 2 opacity uncertainties or photometric errors have biassed the 
results. 

Seven white dwarfs, half of our sample, have 4120 K < T cS < 4480 K and cooling ages 
between 7.3 Gyr and 8.7 Gyr. Their distances are 60-80 pc, and their tangential velocities 
are 40-85 kms -1 . The main-sequence lifetime of these remnants is likely to have been ~2 
Gyr (e.g. Catalan et al. 2008), hence these white dwarfs are most likely to be thick disk 
10-11 Gyr-old objects. Of these seven, we could not constrain the atmospheric composition 
for two, two have mixed atmospheres, two pure hydrogen and one pure helium. We note 
that there are no pure helium atmosphere white dwarfs cooler than 4400 K in the Kilic et 
al. (2010a) sample. The uncertainty in the photometry would allow ULAS J1006+09 to 
be hydrogen-rich, in which case it is slightly cooler, with T eff = 4230 K, and older, with a 
cooling age of 8.4 Gyr. 

Four of the white dwarfs have 4610 K < T e fj < 4810 K and cooling ages between 6.5 Gyr 
and 6.9 Gyr. Their distances are 65-100 pc, and velocities are 60-95 kms" 1 . The remaining 
three white dwarfs have 4970 K < T cff < 5260 K and cooling ages between 4.3 Gyr and 
5.8 Gyr. Their distances are 60-100 pc, and velocities are 70-100 kms -1 . This half of our 
sample may consist of thin disk remnants with unusually high velocities (e.g. Bergeron 2003, 
Tetzlaff et al. 2011), or lower-mass remnants of thick disk or halo late-F or G stars. 

5. Conclusions 

Our search of the 1400 deg 2 of SDSS and UKIDSS LAS (Data Release 6) sky for white 
dwarfs cooler than 5000 K resulted in 17 candidates. One of these was a previously known 
white dwarf, and three proved to be subdwarf or dwarf G- to K-type stars, as determined 
by GMOS spectroscopy obtained at Gemini Observatory. Of the thirteen newly confirmed 
white dwarfs, two show H a, the remaining are featureless. The relatively high number of 
subdwarfs in the sample is due to errors in the initial proper motion derivation, in one case 
due to an error in epoch determination, and in another to the binary nature of the source. 
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We fit the SDSS and LAS ugrizYJHK photometry using model atmospheres which 
can be calculated with arbitrary amounts of hydrogen and helium. As parallaxes are not 
known for these objects, we adopted the canonical white dwarf surface gravity of log g = 8.0, 
although we explore the impact of varying gravity. We have found that in searching for 
objects faint at H, objects are scattered into our sample where the photometry is in error by 
twice the estimated uncertainty. Repeat H- and K-band photometry proved to be necessary 
to fit the white dwarf energy distributions, where H, K > 18 magnitudes. 

Seven of the newly identified white dwarfs have 4120 K < T c g < 4480 K, and cooling 
ages between 7.3 Gyr and 8.7 Gyr; these are likely to be thick disk 10-11 Gyr-old objects. 
The other half of the sample has 4610 K < T e g < 5260 K, and cooling age between 4.3 
Gyr and 6.9 Gyr. These are either thin disk remnants with unusually high velocities (e.g. 
Bergeron 2003, Tetzlaff et al. 2011), or lower-mass remnants of thick disk or halo late-F or 
G stars. Our earlier search of 280 deg 2 of LAS sky (Lodieu et al. 2009) resulted in seven 
spectroscopically confirmed new white dwarfs, with effective temperature 5400 K < T c g < 
6600 K and cooling age 1.8-3.6 Gyr; these are a more extreme example of apparently young 
objects with high tangential velocities. 

We have determined that pairing the SDSS and LAS databases and using proper motion 
and color selections of: H g > 20.5, r < 20.7, 14 < J < 19.6, H < 18.9, 0.9 < g - r < 1.6, 
0.2 < r - i < 0.6, 0.6 < i - J < 1.4, J - H + (y/err(J) 2 + err(H) 2 ) < 0.2 produces 
a reasonably complete sample of 4200 K white dwarfs, which have cooling ages of around 
8 Gyr. When complete, the LAS should provide a sample of about twenty cool, thick disk, 
white dwarfs, and possibly identify white dwarf remnants of the halo. The results presented 
here for various color selections will be a useful guide for identifying cool white dwarfs in 
ongoing near-infrared sky surveys such as UKIDSS and the Visible and Infrared Survey 
Telescope for Astronomy (VISTA; McPherson et al. 2006). This is particularly true for 
white dwarfs showing significant pressure-induced H 2 opacity, due to either extremely low 
temperatures or higher-pressure mixed H-He atmospheres. The near-infrared surveys will be 
useful for studying the complex spectral evolution of white dwarfs, by providing the ratio of 
hydrogen- to helium-rich white dwarfs at different values of T c g. 

We are very grateful to the referee, M. Kilic, for a review that greatly improved the 
paper. 

Some of the data reported here were obtained as part of the United Kingdom Infrared 
Telescope (UKIRT) Service Programme; UKIRT is operated by the Joint Astronomy Centre 
on behalf of the Science and Technology Facilities Council of the U.K.. This paper makes 
extensive use of the UKIRT Infrared Deep Sky Survey (UKIDSS); we are grateful to the 



-14- 



WFCAM instrument team, the UKIRT staff, the UKIDSS team, the CASU data processing 
team, and the WSA group at Edinburgh, all of whom have contributed to the success of the 
survey. 

This paper also makes use of Sloan Digital Sky Survey (SDSS) data. Funding for the 
SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating 
Institutions, the National Science Foundation, the U.S. Department of Energy, the National 
Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck 
Society, and the Higher Education Funding Council for England. The SDSS Web Site is 
: / /www.sdss.org/j The SDSS is managed by the Astrophysical Research Consortium for 
the Participating Institutions. The Participating Institutions are the American Museum of 
Natural History, Astrophysical Institute Potsdam, University of Basel, University of Cam- 
bridge, Case Western Reserve University, University of Chicago, Drexel University, Fermilab, 
the Institute for Advanced Study, the Japan Participation Group, Johns Hopkins University, 
the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and 
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Fig. 1. — Reduced g-band proper motion as a function of g — i. White dwarf cooling 
curves for different tangential velocities are shown; the v tan =25-40 kms -1 curves represent 
the expected location of disk white dwarfs, and the f ta n =150 kms -1 curve the halo white 
dwarfs. The typical location of subdwarfs and stars are also shown (based on Figure 6 of 
Kilic et al. 2010a). Filled circles are confirmed white dwarfs and open circles are G- to K- 
type (sub) dwarfs identified in this work. Grey triangles are the warmer sample of confirmed 
white dwarfs from Lodieu et al. (2009b). Asterisks are cool white dwarfs from Kilic et al. 
(2010a; white dwarfs identified both by the RPM method and spectroscopically by Gates et 
al. (2004) and Harris et al. (2008) are shown). Points circled in red are white dwarfs cooler 
than 4000 K. The regions labelled A and B indicate the search criteria used here, see Figure 
2 and the text. The photometry is taken from DR8 of the SDSS, and the proper motions 
are derived from the DR8 releases of the SDSS and LAS. 
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Fig. 2. — Plots demonstrating the color selection of the white dwarf candidates. Model 
sequences for white dwarfs with pure hydrogen and pure helium atmospheres are shown as 
solid and dashed lines, and a sequence with He/H=100 is shown as a short-long-dash line. 
T ef j/1000 K is indicated along the sequences for the pure hydrogen case. Symbols are as in 
Figure 1. The regions labelled A and B indicate the search criteria used here, as described 
in the text. The optical photometry is taken from SDSS DR8, and the infrared from our 
WFCAM and NIRI observations where available, otherwise from DR8 of the LAS. 




Fig. 3. — Plots of the red to far-red colors of the samples considered here. Symbols, sequences 
and data sources are as in Figures 2 and 3. 
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Fig. 4. — GMOS spectra of the three ULAS G and K stars in our sample, compared to 
similar spectra taken from Le Borgne et al. (2003; smoothed to the resolution of our data). 
The spectra have been normalized to unity at 600 nm, and offset by 0.5 flux units for clarity. 
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Fig. 5. — GMOS spectra of previously unknown white dwarfs in our sample. The spectra 
have been normalized to unity at 600 nm, and offset by 0.5 flux units for clarity. H a 
absorption is seen in the ULAS J1323+12 spectrum, and weakly in the ULAS J2331+15 
spectrum (see also Figures 6 and 7). 
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Fig. 6. — Two of the four white dwarfs in our sample best fit with pure hydrogen atmospheres. 
The error bars in the left panels represent SDSS ugriz and NIRI, WFCAM and/or LAS 
YJHK photometry. The dashed error bars indicate photometric datapoints that have been 
ignored in the fits: u data, as well as g for ULAS J0826— 00, and H and K for ULAS 
J1323+12 (see text). The solid lines in the left panels represent the monochromatic model 
fluxes while the solid dots correspond to the average over the filter bandpass. A surface 
gravity log g = 8.0 is assumed, and the derived T e g is shown in the legends. The right panels 
show the observed spectrum around H a, with the modelled pure-hydrogen atmosphere line 
profiles. ULAS J0826— 00 is too cool to show H a. 
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Fig. 7. — The two other white dwarfs in our sample best fit with pure hydrogen atmospheres. 
Symbols are as in Figure 6. u data have been ignored in the fits, as well as g and K for 
ULAS J1345+15 (see text). ULAS J1345+15 is too cool to show H a (right panels). 
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Fig. 8. — Four white dwarfs in our sample best fit with pure helium atmospheres. Symbols 
are as in Figure 6. u data have been ignored in the fits, as well as H for ULAS J1006+09 
(see text). All objects have observed and modelled featureless spectra (not shown). 
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Fig. 9. — One white dwarf in our sample best fit with either a pure helium atmosphere 
(left panel) or a mixed atmosphere with a high He/H ratio (right panel); the mixed fit is 
preferred. Symbols are as in Figure 6. u data have been ignored in the fits (see text). Both 
the observed and modelled spectrum (not shown) is featureless. 
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Fig. 10. — Two white dwarfs in our sample best fit with mixed H/He atmospheres. Symbols 
are as in Figure 6. u data have been ignored in the fits (see text). Solutions with two 
values of the He/H ratio are possible (see text). Both objects have observed and modelled 
featureless spectra (not shown). 
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Fig. 11. — Two white dwarfs in our sample for which composition could not be constrained, 
the pure hydrogen fit is shown in the left panels, and the pure helium fit in the right panels. 
Symbols are as in Figure 6. u data have been ignored in the fits, as well as g in the pure 
hydrogen fits. Both objects have observed and modelled featureless spectra (not shown). 
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a Uncertainty in proper motion is ~ 14 mas yr 1 , except for ULAS J1323+12 for which the uncertainty is 30 mas yr 1 . 

b Other sources: (1) SDSS-USNO proper motion reported in DR8 of the SDSS; (2) derived here from digital sky images together with the SDSS and LAS 
astrometry; (3) derived here from digital sky images, as well as the USNO-B, SDSS and LAS astrometry. Typical uncertainty is ~ 10 mas yr" 1 . 

c Valuc in brackets is that implied by the alternative proper motion value. 

d All region selections include r AB < 20.7 and 14 < J < 19.6; in addition for region A H g > 20.5, H < 18.9, 0.8 < g - r < 1.6, 0.2 < r - i < 0.6, 
0.6 < i — J < 1.4, J — H < 0.2. Region Ab is the same, except that the source is slightly bluer: g — r = 0.63. Region B is defined by H g > 20.5, 0.2 < g — r < 1.2, 
-0.6 < r-i < +0.6, J-H < -0.1, H-K < -0.1. Region C is defined by H g > 22.5, -0.5 < g-i < 1.0, 0.0 < r-i < 0.6, H < 18.8, J-H < 0.1, H-K < 0. 

e Confirmed as a white dwarf spcctroscopically in this work. 

' Confirmed as a G- to K-typc (sub)dwarf star spectroscopically in this work. 

g Erroncously high initial proper motion determination placed object in sample. 



h This close pair was incorrectly matched, resulting in a high initial proper motion determination which placed the Northern object in sample. 
'No optical spectrum obtained. 

j Previously discovered in SDSS by Kilic ct al. (2010a). 
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Table 2. SDSS DR8 and UKIDSS LAS DR8 Photometry for Candidate White Dwarfs 



Short Name 


u(err) 


3 (err) 




r(err) 




i(err) 


i 




z (err) 


Y (err 


) 


J (err) 


Jf(err 


) 


ii"(err 


) 


ULAS 


J0121-00 a - b 


22 


87(0 30") 


20 


78(0 031 




. ( O^U.UZ f 


1 


o t \yj 


091 
UZ ) 


1 Q 
iy 


1 7(0 0.A1 




OA 1 ! 


J-O.OU^U.UO ) 


lo.OU^U 


is; 


J-O.UU^U. 


9^1 
ZO j 


ULAS 


J0826-00 a - b 


23 


.19(0.44) 


21 


02(0 041 


zu 


f\7(C\ 0.91 
U 1 ^U.UZ J 


1 


7°J fl 


f19l 
UZ ) 


iy 


a\(c\ n^l 


ry.uo^u. 


u ' ) 


lO.UO^U. -LU ) 


lo. ( o^u 


1 71 
.LI) 






ULAS 


J0840+05 a - b 


22 


.37(0.21) 


20 


.61(0.02) 


1 Q 
±y 


Q7(C) 0.91 
-o 1 ^U.UZ J 


1 




■ UZ f 


1 Q 
ry 


-OU^U.UU ) 


in 9n/'n 


UU ) 


j_y .uz^u.uo ) 


in 91 f'n 
ly.z l 


is; 






ULAS 


J1006+09 a 


23 


.26(0.48) 


21 


.11(0.04) 


on 
zu 


CtAiCl n.91 

-U^U.UZ j 


1 Q 

±y. 


*\A((\ 

.O^U 


n.91 

.UZ ) 


1 Q 
±y 


aq(o nsl 


Lo. i oyj. 


UU ) 


1 £ AAff\ f\7\ 
lo.^t^t^U.U I ) 




1 Ql 


1 7 QiVn 


1 ^il 
Li) 


ULAS 


J1142+00 c 


20 


.47(0.04) 


19 


.41(0.01) 


±y 


i nYn nil 


1 Q 

±y. 


ni (i\ 


m 1 

VL) 


1 Q 
iy 


m (c\ n^l 


1 8 A7(l\ 


UU ) 


lo. -LU^U.U i ) 


lo.UU^U 


uy j 


17 Q^/'n 


1 ^1 
LO) 


ULAS 


J1206+03 a - b 


23 


.30(0.52) 


21 


.09(0.04) 


20 


1 t(Cl 0^1 

. J-O KJO J 


19 


65(0 


.03) 


19 


^(C\ OKI 


1 8 QQfO 


08) 


1 & 79('n 1 Vl\ 

-LO. / xi^U. -LU j 




22) 


1 S 4fi('0 


22) 


ULAS 


J1234+06N C 


22 


.74(0.31) 


21 


.02(0.03) 


20 


.14(0.02) 


19 


.81(0 


.02) 


19 


.62(0.06) 


19.11(0. 


05) 


18.81(0.07) 


19.66(0 


26) 


17.95(0. 


13) 


ULAS 


J1234+06S d 


25 


.37(0.65) 


22 


.43(0.10) 


21 


.13(0.05) 


20 


.57(0 


.04) 


22 


.27(0.44) 


19.61(0. 


08) 


19.08(0.09) 


19.28(0 


.18) 






SDSS J1247+06 c 


20 


.95(0.08) 


20 


.03(0.02) 


18 


.68(0.01) 


18 


.39(0 


.01) 


18 


.27(0.02) 


17.76(0. 


02) 


17.54(0.03) 


17.49(0 


04) 


17.54(0. 


10) 


ULAS 


J1320+08 a > b 


21 


.56(0.14) 


20 


.56(0.03) 


20 


.02(0.02) 


19 


.78(0 


.02) 


19 


.70(0.09) 


19.21(0. 


06) 


19.03(0.08) 


19.53(0. 


.16) 






ULAS 


J1323+12 a 


21 


.91(0.17) 


20 


.32(0.02) 


19 


.73(0.02) 


19 


.50(0 


.02) 


19 


.25(0.06) 


18.85(0. 


05) 


18.57(0.06) 


18.64(0 


17) 


18.25(0. 


16) 


ULAS 


J1345+15 a 


23 


.52(0.50) 


20 


.78(0.03) 


19 


.61(0.02) 


19 


.18(0 


.01) 


18 


.95(0.04) 


18.30(0. 


03) 


18.07(0.04) 


17.99(0 


07) 


18.24(0. 


15) 


ULAS 


J1351+12 a - b 


22 


.66(0.29) 


20 


.98(0.03) 


20 


.06(0.02) 


19 


.72(0 


.02) 


19 


.58(0.06) 


19.08(0. 


.07) 


18.79(0.08) 


18.88(0 


.16) 


18.63(0. 


19) 


ULAS 


J1404+13 a - b 


20 


.79(0.07) 


19 


.62(0.01) 


19 


.04(0.01) 


18 


.84(0 


.01) 


18 


.92(0.04) 


18.71(0. 


04) 


18.81(0.07) 


19.01(0 


15) 






ULAS 


J1436+05 a - b 


22 


.05(0.17) 


20 


.29(0.02) 


19 


.38(0.01) 


18 


.87(0 


.01) 


18 


.44(0.03) 


17.85(0. 


03) 


17.62(0.04) 


17.51(0 


06) 


17.47(0. 


.11) 


ULAS 


J1454-01 a - b 


21 


.86(0.12) 


20 


.01(0.02) 


19 


.12(0.01) 


18 


.79(0 


.01) 


18 


.65(0.03) 


18.06(0. 


03) 


17.76(0.04) 


17.68(0. 


.09) 


17.56(0. 


12) 


ULAS 


J2246-00 c 


23 


.28(0.53) 


21 


.35(0.05) 


20 


.66(0.03) 


20 


.28(0 


.04) 


20 


.46(0.16) 


19.54(0. 


.12) 


19.16(0.13) 


19.30(0. 


.35) 






ULAS 


J2331+15 a - b 


22 


.92(0.33) 


20 


.51(0.02) 


19 


.77(0.02) 


19 


.44(0 


.02) 


19 


.25(0.05) 


18.78(0. 


05) 


18.50(0.08) 


18.72(0. 


.25) 







Note. — SDSS ugriz magnitudes are on the AB system (Fukugita et al. 1996). LAS YJHK arc on the Vega system (Hewett et al. 2006). 
a Confirmcd as a white dwarf spectroscopically in this work. 
b Improved H and K photometry is given in Table 4. 

c Confirmed as a G- to K-type (sub)dwarf star spectroscopically in this work. 
d No spectrum obtained. 

Previously discovered in SDSS by Kilic et al. (2010a). 
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Table 3. GMOS Observation Log 



Short Name 


SDSS 


r 


Total Exp. 


Date 


Program 






AB 




seconds 


YYYYMMDD 




ULAS 


J0121-00 a 


19. 


,70 


5100 


20100710 


GS-2010A-Q-58 


ULAS 


J0826-00 a 


20. 


,07 


8400 


20100119 


GS-2010A-Q-58 


ULAS 


J0840+05 a 


19. 


.83 


3600 


20081209 


GS-2008B-Q-35 


ULAS 


J1006+09 a 


20. 


,04 


6400 


20090122 


GN-2008B-Q-111 


ULAS 


J1142+00 b 


19. 


.11 


1800 


20100119 


GS-2010A-Q-58 


ULAS 


J1206+03 a 


20. 


,13 


5040 


20090120 


GS-2008B-Q-35 


ULAS 


J1234+06N b ' c 


19. 


,76 


5760 


20090122 


GS-2008B-Q-35 


ULAS 


J1320+08 a 


20. 


.02 


6400 


20090126, 20090127 


GN-2008B-Q-111 


ULAS 


J1323+12 a 


19 


.77 


4800 


20090122, 20090125 


GN-2008B-Q-111 


ULAS 


J1345+15 a 


19. 


.61 


4400 


20100415 


GS-2010A-Q-58 


ULAS 


J1351+12 a 


20. 


,07 


8400 


20100418, 20100607 


GS-2010A-Q-58 


ULAS 


J1404+13 a 


19 


,04 


1800 


20100508 


GS-2010A-Q-58 


ULAS 


J1436+05 a 


19. 


.40 


3200 


20100607, 20100713 


GS-2010A-Q-58 


ULAS 


J1454-01 a 


19. 


,12 


1800 


20100418 


GN-2008B-Q-111 


ULAS 


J2246-00 b 


20. 


.64 


11400 


20080731 


GS-2008B-Q-35 


ULAS 


J2331+15 a 


19 


,77 


4800 


20081221 


GN-2008B-Q-111 



a Confirmed as a white dwarf spcctroscopically in this work. 

b Confirmed as a G- to K-type (sub)dwarf star spectroscopically in this work. 

c Close pair of objects, spectroscopy is for northern source. 
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Table 4. WFCAM and NIRI Photometry for LAS White Dwarfs 



Short Name 


if (err) 


K(eri) 


Exp. H, K 


Date 


Program 










minutes 


YYYYMMDD 




ULAS 


J0121-00 


18.13(0.02) 


18.05(0.03) 


16, 16 


20100814,20100815 


USERV 1876 


ULAS 


J0826-00 


18.59(0.02) 


18.58(0.02) 


9, 9 


20110126 


GN-2011A-Q-59 


ULAS 


J0840+05 


19.21(0.03) 


19.89(0.06) 


13.5, 18 


20110125 


GN-2011A-Q-59 


ULAS 


J1206+03 


18.46(0.03) 


18.33(0.06) 


16, 16 


20100613 


USERV 1876 


ULAS 


J1320+08 


19.08(0.03) 




22.5 


20110126 


GN-2011A-Q-59 


ULAS 


J1351+12 


18.57(0.03) 


18.48(0.05) 


16, 16 


20100613 


USERV 1876 


ULAS 


J1404+13 


19.18(0.05) 


19.72(0.16) 


16, 16 


20100613 


USERV 1876 


ULAS 


J1436+05 


17.66(0.03) 


17.49(0.03) 


1.25, 1.25 


20110105 


GN-2011A-Q-59 


ULAS 


J1454-01 


17.64(0.03) 


17.49(0.03) 


1.25, 1.25 


20110209 


GN-2011A-Q-59 


ULAS 


J2331+15 


18.28(0.03) 


18.26(0.04) 


16, 16 


20100807, 20100812 


USERV 1876 



Note. - 



HK are on the Vega system (Tokunaga & Vacca 2005, Hewett et al. 2006). 
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Table 5. Derived Properties of the LAS White Dwarfs 



Short Name 


Spectral 


Atmospheric 


T off a 


Cooling b 


Distance 


«tan d 




Type 


Composition 


K 


Age, Gyr 


pc 


km s _1 


TTT A Q TH1 01 flfl 




unconstrained 


A A OA I QOfl 


70 + 2.0 

' - 3.7 


Q1 11 A 


qo_Li n rci 1 
oo±lU [olj 


ULAS J0826-00 


DC 


H 


4170±130 




76±10 


33±9 [46] 


ULAS J0840+05 


DC 


mixed 


4350±110 




82±15 


44±9 [50] 


ULAS J1006+09 


DC 


He 


4260±100 


7 7+O.4 
' • ' — 3 2 


73±17 


51±12 [50] 


ULAS J1206+03 


DC 


He 


4610±80 


6-9±g:i 


94±20 


73±15 [71] 


SDSS J1247+06 


DQpcc 


mixed 


5120 


5.6 


60 


110 


ULAS J1320+08 


DC 


high-He mix 


4800±130 


°-° -3.0 


102±18 


95±19 [101] 


ULAS J1323+12 


DA 


H 


5260±90 


4.311-9 


103±21 


163±30 [102] f 


ULAS J1345+15 


DC 


H 


4120±150 


8-7 1 1 ^ 


58±8 


72±10 [71] 


ULAS J1351+12 


DC 


He 


4810±70 


6.5t° 3 9 2 


103±21 


59±12 [51] 


ULAS J1404+13 


DC 


mixed 


4380±100 


7 r +0.4 
' '° - 2.9 


59±10 


43±8 [40] 


ULAS J1436+05 s 


DC 


unconstrained 


4340±170 


7.8t 2 3 .7 


57±12 


85±15 [81] 


ULAS J1454-01 


DC 


He 


4800±60 


6-5 i° 3 9 2 


65±15 


65±15 [83] 


ULAS J2331+15 


DA 


H 


4970±80 


5-8t 3 3 i 


95±27 


71±20 [67] 



a The uncertainty in T e g is due to both photometric scatter and an allowed range in gravity of 
7.7 < \ogg < 8.3. 

b The cooling age is derived from the composition and temperature using the cooling models of 
Fontaine, Brassard & Bergeron (2001). The uncertainty is due to the range in gravity of 7.7 < log g < 
8.3. 

c The distance is estimated from the modelled and observed magnitudes (Holberg & Bergeron 2006). 
The uncertainty is due to the range in gravity of 7.7 < log g < 8.3. 

d The tangential velocity is calculated from distance and the proper motion given in Table 1. Values 
in brackets are those implied by the alternative proper motion given in Table 1. 

e Previously discovered in SDSS by Kilic et al. (2010a), tabulated properties are from that work. 

f The alternative value for v tan is preferred, given the large uncertainty in the SDSS-LAS proper 
motion caused by the short timeline. 

s The fit to this white dwarf is poor (see Figure 11). 



